Perturbative Spectrum of Trapped Weakly Interacting Bosons in 

Two Dimensions 



Velimir Bardek, Larisa Jonke, and Stjepan Meljanac * 
Theoretical Physics Division, 
Rudjer Boskovic Institute, P.O. Box 180, 
HR-10002 Zagreb, CROATIA 

Abstract 

We study a trapped Bose-Einstein condensate under rotation in the limit 
of weak, translational and rotational invariant two-particle interactions. We 
use the perturbation-theory approach (the large-iV expansion) to calculate 
the ground-state energy and the excitation spectrum in the asymptotic limit 
where the total number of particles N goes to infinity while keeping the total 
angular momentum L finite. Calculating the probabilities of different config- 
urations of angular momentum in the exact eigenstates gives us a clear view 
of the physical content of excitations. We briefly discuss the case of repulsive 
contact interaction. 
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The study of low-lying excitations of the weakly interacting, trapped Bose-Einstein con- 
densed gas under rotation is of considerable experimental and theoretical interest @. 
Theoretical studies have focused on the Thomas- Fermi limit of strong interactions ||, as 
well as on the limit of weak interactions [§-0], which we consider in this paper. Wilkin et al. 
Q studied the case of attractive interaction, and Mottelson and Kavulakis et al. devel- 
oped a theory for repulsive interactions. They compared the mean-field approach and exact 
numerical results obtained by diagonalization in a subspace of degenerate states 0. Bertsch 
and Papenbrock |7j performed numerical diagonalization for small systems and showed that 
the interaction energy of the lowest-energy states decreases linearly with angular momentum 
L. Nakajima and Ueda || found through an extensive numerical study, in the limit where 
the angular momentum per particle is much smaller than one, that low-lying excitation en- 
ergies, measured from the energy of the lowest state are given by 0.795 n(n — 1), where n is 
the number of octupole excitations. Recently, Kavoulakis et al. || rederived these results 
analytically with use of the diagrammatic perturbation-theory approach in the asymptotic 
limit iV — > oo. In this paper we present a systematic method for calculating the excitation 
spectrum for the weak, translationally and rotationally symmetric interaction in the asymp- 
totic limit, where the total number of particles N goes to infinity, while keeping the total 
angular momentum L finite. We also discuss the probabilities of different configurations of 
the angular momentum in the exact eigenstates. 

Our starting point is the two-dimensional Hamiltonian H = H + V, where 

is the one-particle part, including the kinetic energy of the particles, and the potential energy 
due to the trapping potential. The trapping potential is approximated by a two-dimensional, 
isotropic harmonic oscillator with the frequency set to one. The system is in the ground state 
for the motion in the direction of the axis of rotation. The two-body interaction between 
the particles is given by 

^ = f>(| rj -r,|), (2) 

i<j 

where an arbitrary potential v possesses translational and rotational symmetries. We also 
assume that the interaction v is weak. This allows us to work within the subspace of single- 
particle states with no radial excitations 

Vv l (z) = (7m!)- 1 /Vexp(-i|«| 2 ) ) (3) 

where z = x + iy and n is the angular momentum quantum number. The energy lev- 
els and the corresponding wave functions are found by diagonalizing the interaction V 
in this Hilbert space. Basis functions for the many-body problem are ip(zi, z 2 , . . . , z^) = 
(p(zi, z 2 , ■ ■ ■ , zn) n£Li exp(— ||zi| 2 ), where ip is a homogeneous polynomial of degree L. For 
simplicity, we omit the exponentials from the wave functions. Suitable basis functions for 
such polynomials are given by 
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1 N 

b x ( Zi , Z2 ,..., ZN ) = — E'4 X 4 2 ■ • • (4) 

1' 2- p- i 1) i 2 ,...,i g =l 

where the set {Ai, A2, • • • , \} denotes any partition of L such that J2f=i K — L and Ai > 
\ 2 > ■ ■ ■ > X q > for q < N Jl(| . The prime denotes the sum over mutually different 
indices i±, i 2 , ■ ■ ■ , i q , while the numbers ui, u 2 , . . . , v v denote the frequencies of appearance 
of equal A's. Note that the number of distinct monomial terms z^z^ ■ ■ ■ z i q in B\ is given 

by d x = N(N - 1) • • • (N - q + 1)/(i/i!i/ 2 ! ■ • • v p \), where V\ + P 2 -\ h v p = q. Owing to 

the translational and rotational symmetries of the two-particle interaction v{\xi — r 3 -|), for 
non-negative integers n and m we have 

v(\ Zl -z 2 \) (z x + z 2 ) n ( Zl - z 2 ) 2m P(z 3 , z 4 , ... , z N ) 
= c 2m {zi + z 2 ) n (zx - z 2 ) 2m P{z 3 , z 4 , ■ ■ • , z N ), (5) 



where P denotes an arbitrary polynomial. The coefficient c n is given by 

joo dr r 2n+l v ^ exp (_ r 2/ 2 ) 

° n ~ f™drr 2n + 1 exp(-r 2 /2) ' 



(6) 



and represents the interaction energy v(r) of the relative motion of two bosons in the single- 
particle state r n exp(— r 2 /2) with the angular momentum n. To proceed, let us now define 
symmetric functions of two variables 

bij(zu z 2 ) = ^{z{z J 2 + z{z l 2 ), i> j. (7) 

The action of the potential u(|r* — r 3 -|) on bij is given by 

[f] 

v(\z! - z 2 \)b i j(z 1 , z 2 ) =^2ai-b ki (z 1 ,z 2 ), (8) 

1=0 



where i+j = k + l = n. This restriction is a consequence of the conservation of total angular 

1. 1 



momentum for a rotationally symmetric potential. Also, the coefficients oM 



9 _ A 12 

kl _ z "kl „ c2p cl 

ij 



n C 2pSi P jS l n -2p,2pi (9) 

p=0 



where 



* = &-KN.)' (10) 

r Tt j 

satisfy the summation rule J2i=o a ij — °o, as a consequence of translational symmetry. The 

coefficients a k j in fact represent the two-body matrix element Vijki (see Ref. JO]]) of the 
interaction potential V. Next, by using Eqs.([|) and (Q) we obtain 

fN\ N 

V B r 1 l {z 1 ,z 2 ,...,z N ) = c \^ 2 jB^(z 1 ,z 2 ,...,z N ) for Si = XX ( n ) 
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As a result, B™(zx, z 2 , . . . , zjv) is an exact eigenstate and the corresponding eigenvalue is 
Co (^j ■ Furthermore, owing to translational invariance the action of the potential V on the 
product B™B\ reduces to 

VB?Bx = B?VB x (12) 

for any n and partition A. Specially, if A is an eigenstate with energy E, then B™A is also 
an eigenstate with the same energy. 

Generally, for any partition A we find 

VB X = Y,< B ^ ( 13 ) 

where /1 is the partition obtained by substituting a pair of numbers {k,l} for a {i,j} in 
partition A, such that i > j , k > I, and i+j = k+l, for all distinct pairs {i, j} and all allowed 
{k, I}. Note that any partition A = {Ai, ■ • • , AJ can be written as A = n °l ni 2 n2 ■■■l n ' 
with Y,i n iW — N and J2i i n iW — L. In the second quantized approach, the number 7ij(A) 
can be interpreted as the number of particles with the angular momentum % in the partition 
A. The diagonal coefficient is 

«l = E4^.w[«iW-y, (14) 

where the sum goes over all distinct pairs {i,j}, % > j > contained in partition A. The 
sum contains ( K ^ + K2 terms, where K\ is the number of n^'s greater than zero , and K2 
is the number of rij's greater than one, in partition A (i > 0). The nondiagonal coefficient 
can be expressed as 

< = <^^M*^)-M, (is) 

where {i,j} ({k, I}) are contained in partition A (/i), respectively. The general matrix 
element has the form = ^-Sx^N 2 + fix^N + 7a m . Note that this matrix is not Hermitian 
since our initial basis {Bx} is orthogonal but not orthonormal, i.e., (B X \B X ) = dxHiXil. 
Since the interaction is Hermitian, changing the basis to orthonormal would render the 
matrix {a^} Hermitian. The matrix {a^} has dimension V(L), which is the number of 
partitions of L. It had been shown [|l2| , |l4f that the eigenvalue problem can be reduced to 
V(L) — ViL — 1) — 1 dimensions and recursively solved for the general interaction up to 
L = 5. For L = 6, the problem reduces to the diagonalization of the 3x3 matrix, which 
can be accomplished using the 1/N expansion. Motivated by this approach, in this paper 
we propose a similar strategy. In the limit where the angular momentum L is much smaller 
than the number of particles N we use perturbation theory in the large-iV expansion to 
calculate the interaction energies and derive analytical results. In the zeroth order, the 
standard perturbation-theory approach gives Ax = Bx for eigenstates and = for the 
corresponding energy. The eigenenergy with the first-order corrections is 
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The above expression is applicable if the condition a^a^ « (a* — a^) 2 is satisfied for all 
partitions // ^ A. It can be easily checked using relations flUD and ([TJ) that a^a^ ~ A" and 
(a^J — a^) 2 ~ A^ 2 . One finds that the dominant contributions are those with j or I equal to 
zero, in Eq. (|T5|) , and they produce corrections to the energy of order iV°. 

We can label the exact interaction energies and eigenstates as E\ and A\, respectively, 
such that in the limit iV — ► oo and finite L, the energy E\ goes to E^ and A\ goes to B\. 
For an exact eigenstate A\(N, L), the state B™A\(N, L) is an exact eigenstate Ay(N, L + n), 
where A' = o' no ~ n )l <ni+ "'* ) 2 n2 • • •. According to translational invariance, we obtain the exact 
identity for eigenenergies 

Hence, Ao"oi n i2' i 2... — B"MQ( no + ni )2n 2 ... and -£7o n oi n i2 n 2... = -^o( n o+ n i)2"2---- The part l ni in the 
partition A denotes n\ unit angular momenta which can be realized only as the angular 
momenta due to the center-of-mass motion. Therefore we consider only the eigenstates 
with partition A = n °2 n2 3 ri3 • • • l ni ■ ■ -, i.e., the states involving quadrupoles, octupoles, and 
higher I poles. In this case, we have 



Er=co[ U ; V +$>!fo(A)no(A)+ £ a% — 7 p^n i (\)[n j (\) - 8 itj ]. (18) 

V A ) i>2 i>j>2 / 



For special partition A = O^ -1 -*/ 1 , we obtain the excitation energy for a general weak 
interaction q = E^' — Eq ^ = A^a^ — a^). In the case of contact interaction it reduces to 
the ei = —c N(l — 2~( / ~ 1 )) ||. Now, we include corrections 

4 1) = 4 0) + £ c^(A)h(A)-5 M ]K +J (A) + l] 

i>j>2 

- J2 CijUi+jiX^miX) + l + 5 iij ][n j (\) + 1], (19) 

i>j>l 



where 



i+j,0 ij 2-<5. 
a ij a i+j,0 2 



c ..- ZM *+J' u 2 (20) 

13 n i0 i fJO _ i+j,0 _ 00 • ^ 
a i0 ' a j0 a i+j,0 "00 

In the case of repulsive delta interaction, Eqs.flT8|) and (|T9|) simplify significantly, because 
all coefficients c n are zero for n ^ 0, so afj = 2^~ n ^(2 — 5k,i)(fyco, n = i + j = k + l. From 
Eq.(|T8"D we easily find the lowest-order energy E 



(0) 
A 



L + 2 _ (i 1 



4 ^ V4 2 l 

i>A 



(21) 



One can calculate the first-order correction for arbitrary partition, so, for example, the 
eigenenergy for the partition A = n °2 n2 3 n3 4 n4 is 

E? = c ^-j(2L + A-n 4 ) + ^n 3 (n 3 - 1) 

/81 27 99 93 M , n/ ; , rX , s 

+ n 4 — n 2 H n 3 H n 4 H H+0(l/iVV (22) 

\52 2 41 194 388/ J V 7 ; V ; 
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This result is in complete agreement with the results obtained in Refs. P,|r2|. We see that 
in the special case of n; = for I > 4, the zeroth-order energy 



E 



(0) 



c 



N 2 N(L + 2) 
~2 4 



is degenerated, but the corrections ^713(713 — 1) remove this degeneracy if n 2 > 2. Hence, for 
the repulsive delta interaction the ground state is unique and defined by N = n + n 2 , L — 



2n 2 or N = no + n 2 + 1, L 



2n 2 + 3, depending on L being even or odd. Therefore, our 

in the limit of large N and finite 



analysis confirms the conjecture of Smith and Wilkin JTTJ 
L. 

The exact eigenstates for L < 5 do not depend on details of interaction, and can be 
expanded in the standard basis B\. For example, in the L = 2 case we have B\ = B 2 + 2Bu 



and A? 



N 



B 



11 



N- 



2N 



-B 2 and in the L = 3 case we have 



B\ 
B X A 2 

A, 



B3 + 3B 2 i 
N — 1 



65m 
N- 



~ Bi ~ '~2PT B21 + N Bim 



2N 

(N -l)(N-2) 
37V2 



N-l 



B 



12 



21 



(N-l)(N-2) 



B x 
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It is interesting to consider the probability of the configuration B^ in the exact eigenstate 
Ax, since it gives us a physical picture of the excitations. One can easily calculate the 
probability using the formula 



w^(A } 



(B, 


A x ) 2 


(B, 


B,)(A X 





(23) 



For the L = 3 case, the probabilities are given in Table 1. We see that the probability of 
the configuration B x l in any exact state other than A\l tends to zero, in the limit when 
iV — > 00 and L finite. Of course, this is not a surprise, as we labeled the exact states to 
obey the condition w^{A x ) — > in the large- N limit. 

The eigenstates cl = Y^{ z h — B\/N) • • • (zi L — B\/N), with L < N are common to all 
interactions and these are the ground states for the repulsive delta interaction ]T3,rj|. There 
is a simple relation between these states and the exact eigenstates A\ we have discused up 
to now. Generally, = A 2 l/2 for even L and = A 32 (l-3)/2 for odd L. For L « N e L is 
dominated by {2 L / 2 } or {32^ L ~ 3 ^ 2 } configurations, depending on L being even or odd. 

Of special interest are "vortex" states ei=N- The probability that every particle carries 
a unit of angular momentum is easily calculated for low N and is given by wn(e 2 ) = 1/2, 
wri(e 3 ) = 4/9, 1014(6:4) = 15/32, wis(e 5 ) = 296/625. It seems that iujjv (ejv) ~ 1/2, an d this 
is in contrast with the naive expectation that all particles contribute with a unit of angular 
momentum in the vortex state [11|. Note that in this case our perturbative approach is 
not valid, hence the probabilities of configurations {2 L / 2 } and {32^ L_3 ^ 2 } are small. 

In conclusion, we studied a trapped Bose-Einstein condensate under rotation in the 
limit of weak, translational and rotational invariant two-particle interactions. We have used 
the perturbation-theory approach to calculate the ground-state energy and the excitation 
spectrum in the asymptotic limit where the total number of particles iV goes to infinity while 
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keeping the total angular momentum L finite. Calculating the probabilities of configurations 
in the exact eigenstates A\ gives us a clear view on the physical content of excitations 

A\. In addition, we have briefly discussed the case of repulsive delta interaction. 
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TABLES 



-Bill B2I B; 



Bf 
BiA 2 
A 3 



NJ V NJ N V NJ N 2 

fi-AV »fi-i 

N\NJ \ N J N \ N 

— -(1--) fi-iVi-- 

a 2 iv v a; 1 aH iv 



TABLE I. The probabilities of configurations -B^ in exact states A^ for the L = 3 case. 
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